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Abstract: We study the shape dependence of entanglement entropy (EE) by deforming 
symmetric entangling surfaces. We show that entangling surfaces with a rotational or 
translational symmetry extremize (locally) the EE with respect to shape deformations that 
break some of the symmetry (i.e. the 1st order correction vanishes). This result applies 
to EE and Renyi entropy for any QFT in any dimension. Using Solodukhin’s formula in 
4d and holography in any d, we calculate the 2nd order correction to the universal EE for 
CETs and simple symmetric entangling surfaces. In all cases we find that the 2nd order 
correction is positive, and thus the corresponding symmetric entangling surface is a local 
minimum. Some of the results are extended to free massive fields and to 4d Renyi entropy. 
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1 Introduction 


Entanglement entropy is a measure of the quantum correlations of a system. It has a 
very wide range of applications from condensed matter physics to quantum field theory 
[1-33] and AdS/CFT [35-62]. The EE in QFT is generally hard to calculate, and most 
computations have been done in simple setups, such as: free fields, CETs, and symmetric 
entangling surfaces (e.g. spheres and planes). Thus there is a need to obtain analytical 
results for interacting theories, non-CETs, and less symmetrical entangling surfaces. This 
work aims to make a step in this direction by studying the shape dependence of EE. 
Previous works on the shape dependence of EE include [63-81]. 

The divergent structure of entanglement entropy for a CFT in d-dimensions is: 

, [ cif+ , d = odd 

+ + ^ d = euenf 


where R is the scale of the entangling region, and 5 is the UV cutoff. The leading divergence 
is the area law, and all of the power law divergences are non-universal. We will be interested 
in the universal term which in d = even is the coefficient of the log divergence and 

in d = odd it is the finite term. 
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Figure 1. Illustration of a perturbed circle r(0) = 1 + e a„ cos(n(/)). Left: A Perturbation 
without a zero mode oq = 0. Right: A Perturbation with a zero mode ag = 0.3. 


Consider a QFT parametrized by coordinates {t,yi,r), where i = 1... ,d — 2. For 
instance take r to be the radial coordinate in spherical coordinates, and yi to be angles 
parameterizing the entangling surface. We will always work in a constant time slice t = 0. 
Consider a codimension-2 entangling surface defined by: 


r{yi) = ro{yi) ( 1 . 2 ) 

where rQ(yi) is some given function of y^. Thus we choose r to be the dependent coordinate, 
and yi as the independent coordinates. We denote the entanglement entropy corresponding 
to the entangling surface ro{yi) as Sq. Now we slightly perturb the entangling surface: 

riVi) = roiyi) + ef{yi) (1.3) 

where e is a small parameter, and f{yi) is some arbitrary perturbation function^. 

The entanglement entropy will change as a result of the perturbation of the entangling 
surface, and it can generally be written as an expansion in e: 

S = So + Sie + S2€^ + ... (1.4) 


The above procedure was carried out in [63, 64] for the case of a perturbed sphere for 
a CFT in d dimensions. They start with a sphere entangling surface in a flat space-time 
background, and perturb the sphere as follows: 


r{ild-2) 




(1.5) 


where R is the sphere radius, the a’s are constants, and the T’s are (real) hyper-spherical 
harmonics. Then they calculate the resulting change in the universal term of the EE. They 
find that the the change in the universal EE vanishes at 1^* order, namely 

^As an example, a perturbed circle in d = 3 is shown in Fig. 1. In this case, (1.3) is given by: r{<p) = 
i?[l + eX)n cos{n<j))], where we Fourier expanded the perturbation /. R is the radius of the circle, (j) is 
the angle in polar coordinates, and a„ are the Fourier coefficients. 
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Q(univ) 


= 0 


( 1 . 6 ) 


Thus, for a CFT the sphere is a local^ extremum with respect to perturbations of 
the entangling surface. Additionally, for holographic CFTs [64] uses the Ryu-Takayanagi 
formula [35, 36] (more specifically, its generalization to higher derivative gravity [83-86]) 
to calculate the 2”“^ order correction 


g{univ) _ 


TT 2 (d — 1) 

2'='-2F(d + 2)F(d) ^ 


2 


■rud- 




k=l 


d = odd 
d = even 


(1.7) 


where Ct is the (positive) central charge appearing in the 2-point function: (TT) ~ 

A priori, could have depended on the three parameters Ct,t 2 -, and t 4 of the 3-point 

function (TTT), but it turns out that it depends just on Ct- 

gf^niv) gigg^j-iy positive, implying that the sphere is a local minimum for holographic 
CFTs. (1.7) was compared to obtained from Solodukhin’s formula in 4d CFTs, and 

precise agreement was found. 

In this work we will generalize the above results of [64] to less symmetric entangling 
surfaces, and (in some cases) to non-CFTs. We will show that entangling surfaces with a 
rotational or translational symmetry in some direction (see Fig. 2), extremize the universal 
EE with respect to shape deformations that break some of the symmetry^, i.e: 

giuniv) ^ 


The proof of this result will be purely geometrical, and hence will apply to any QFT (and 
also to Renyi entropy). A simple corollary is that for d = even, (1.8) is true also for 
multiply connected entangling surfaces^. 

Additionally, we will calculate the 2'^'^ order correction (whose sign determines if 

the local extremum is a minimum or a maximum) for some simple entangling surfaces. We 
perform this calculation using holography (the Ryu-Takayanagi formula), and also by using 
Solodukhin’s formula (for 4d CETs). In all of the the examples of symmetric entangling 
surfaces that we have checked, the 2”'’* order correction is positive, and this corresponds 
to a local minimum. We conjecture this to hold more generally to symmetric entangling 
surfaces. We also comment on results for free massive fields, and 4d Renyi entropy. 

the topology of the entangling surface is allowed to change, then the EE can become unbounded from 
below, as shown in 4d in [82]. 

®It should be emphasized that the family of shape perturbations that extremize the EE, are only those 
which break some of the symmetry of the original entangling surface. For example, for a 3d circle the 
perturbation in Fig. 1-Left is allowed, whereas the perturbation in Fig. 1-Right is not allowed because it 
has a component in the radial direction (a zero mode of the Fourier expansion). Such perturbations cannot 
deform a sphere into a larger/smaller sphere. More generally, such perturbations cannot deform a surface 
of revolution into a different surface of revolution (containing the same number of symmetries). 

^For d = even the universal term of the EE is a log divergence which is determined locally by the shape 
of the entangling surface. Thus for a multiply connected entangling surface the log term is a superposition 
of the contribution from each separate piece of the entangling surface. 
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Figure 2. Left: An example of a d = 4 surface of revolution entangling surface with rotational 
symmetry. Right: A d = 4 waveguide entangling surface with translational symmetry. Such 
symmetric entangling surfaces can obviously be generalized to higher dimensions. 

2 The First Order Correction: Stationarity 

Consider a QFT parametrized by coordinates {4>, yi, r), where i = 1... ,d — 3. Assume (p 
to be a symmetry direction of the entangling surface. Now perturb the entangling surface 
(see (1.3)) with a single Fourier mode^ : 


r{4>, yi) = ro(yi) + eAnan{yi) cos(n(/>) , n / 0 (2.4) 

where ro(?/i) doesn’t depend on the symmetry direction cj), the an{yi) are functions of yi, 
and An are constants. The resulting EE can be expanded as: 


5 = 5o + e5i + e^S2 + ... (2.5) 

Now lets consider the same perturbation but with negative sign, i.e. e —>■ —e: 

f((/), yi) = ro{yi) - eAnUniyi) cos{n(p) , n / 0 (2.6) 

®A general perturbation (without a zero mode) can be written as: 

r{(p, yi) = roivi) + [A„a„{yi) cos{n<p) + B„bn{yi) sin(n((>)^ (2.1) 

n^O 

Then the EE can be expanded: 

S' = 5o + eS'i + e^Sa + ... (2.2) 

At linear order in e the modes don’t mix, and their contributions to Si add up linearly: 

Si = ^ (A„[/i(n) + BnU2{n)') (2.3) 


where Ui{n) and U 2 (n) are some functions of n. Because the modes don’t mix at linear order, we can 
compute the contribution of a single mode and then sum over all modes. 
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The resulting EE will be (just flipping the sign of e in (2.5): 


S = So-eSi + e^S 2 + ... 


(2.7) 


But the two perturbations (2.6) and (2.4) describe precisely the same entangling sur¬ 
face, only rotated. This can be seen by performing ^ on (2.6), which gives (2.4). 

Since the two entangling surfaces are the same, they have same EE and therefore from 
(2.5) and (2.7) we have: 


( 2 . 8 ) 


S = S 


Si = Q 


and we proved what we wanted. 

An important point to make is that the proof above used only the rotation symmetry, 
and did not use any specific property of entanglement entropy or of the QET. Therefore 


stationarity will hold for any quantity which is a function of the surface (e.g. Renyi 


entropy), and for any QET. 

3 The Second Order Correction in Holography 

Consider a boundary QET parametrized by coordinates (t, r), z = 1..., d — 2. z is the 

holographic coordinate, and we set t = 0. The holographic EE according to the Ryu- 
Takayanagi formula is: 



(3.1) 


where £ = . The corresponding equation of motion for the bulk surface is: 


_ d dC _ d dC 
dr dz d{dzr) dyi d{dy^r) 


(3.2) 


where there is a summation convention on yi. 
Consider an entangling surface defined by: 


r{yi) = roiVi) 


(3.3) 


Now perturb the entangling surface: 


r{yi) = ro{yi) + ef{yi) 


(3.4) 


where e is a small parameter. Then the bulk surface will also get perturbed®: 
r{z, yi) = roiz, yi) + eri{z, yi) + e^r 2 {z, yQ + ... 


(3.5) 


®We use the same letter r for the entangling surface and the bulk surface. They can simply be distin¬ 
guished by the fact that for the bulk surface there is a « dependence. 
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For the rest of this section we assume d = even dimensions, and we will want to 
compute the universal log term. We can then write the resulting perturbed C and S as: 

C,Q + C,\e + C-2^ + • • • 

(3.6) 

We can derive a formula for the 2 ’^'^ order correction S^- Assuming that the entangling 
surface has a symmetry in all directions (e.g. the cylinder surface in the next section), 
terms containing r 2 will fall and we get: 


5 = 5o + Sie + 526^ + ... = 


d^ ‘^UidzC = 


d'^ “^Vidz 



d-2 


Vidz 


ri 


d d dC 
de dm d{dyT 


+ dyji 


d dC 

ded{dy,r) 


+ / 


d dC 


e=0 


IBM 


de d{dzr) 



^max 1 

0 

II 

z=S ) 


(3. 


3.1 Cylinder Entangling Surface 

Now let us assume that the entangling surface is a cylinder x in flat space- 

time. We denote yj as cartesian directions along (where —L < yj < L) and Vtp 

as directions along the sphere (of radius R). The holographic EE written in these 
cylindrical coordinates is (see also (A. 20 )): 

S = [ dzd'^~‘^~^yjdQp C{z,r,yj,Qp) (3.8) 


where: 


Cyz^r 


1 F 2 {z){dzrY + F^{z) {dy^rf + -^{dn^rY 


(3.9) 


We are not committing yet to a particular metric in the bulk, therefore Fi{z), ^ 2 ( 2 ), ^ 3 ( 2 ;) 


R. 


are arbitrary functions of z. Eor an AdS metric we have: F 2 {z) = F 3 {z) = 1, and Fi = 


We Eourier expand the perturbation / in (3.4): 


d—2—p 

f{yj,np)= ^ \^a{nj,ip}YiY^p) cos(nj7/i) 


(3.10) 


{TLj ^Ip } 2 1 

where Yi^{idp) are (real) hyperspherical harmonics, the nj are integers, the a^nj,ip} are the 
coefficients of the Fourier expansion, and for conciseness we defined {nj} = ni ,..., nd- 2 -p 
and {Ip} = I, mi ..., mp-i. 

The bulk surface perturbation ri of (3.5) can also be Fourier expanded (see (A.25)): 


ri{z,yj,%)= ^ 

{Uj ,/p} 


( 1 ) 

.Ip}^ {uj ,lp} 


d—2—p 

n cos(niyj) 


2=1 


(3.11) 
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where ^ |(z) are functions of z, which must obey the boundary condition: ^ |(z = 

Now we can derive an expression for the boundary term (the last term on the RHS of 
(3.7)): 


„(i) 


bound. _ L^-^-^F,{z)F2{zy, ^ 


v.P -1 


S: 


{Uj ,/p} 


2z'^-i[l + F2(9^ro)2]V2 ^ 

{rij ,/p| 


f) 

^ o „ /^(i) ')2 

l + F2(5zro)^ ^ 


=5 

(3.12) 


where L is the length of the directions, and we take the limit L —?■ oo. This formula 

is a boundary term which is evaluated at the boundary z = 5. We will want to extract the 
universal log divergence from this formula in d = even dimensions. 

In principle, rQ{z) and ^ ^{z) can be obtained by solving the EOM’s for the bulk 
minimal surface (A.23) and (A.26). These solutions will generally have Fefferman-Graham 
type expansions of the form ([19, 43, 87, 88]): 

roiz) = qo + q 2 Z^ + ... + + qdz'^\og{z) + ... (3.13) 


= ^0 + U 2 Z^ + . . . + UdZ^^ + UdZ'^ log(z) + 
The terms in (3.12) can then be written as: 


(3.14) 


„(i) ^ ..( 1 ) pd^ro{r\^J, 


[1 + F2{dzrQy]z‘ 


d-l 


= dqououd log z 


{nj,lp}^ 


yd—1 


log 


= dpulqdlogz 


(3.15) 


where we extracted the log terms. Using (3.12), and also the fact that the background 

d 2 d-\-2 d 

is asymptotically AdS: Ti(0) = = —— r(d+ 2 )^^^^7 ^2(0) = 1, we get^ 


abound. I 

2 I log 


, . d — 2 d-\-2 / y/ \ 

(_!)—^—r(|) 
(d + i)r(d-i) 


CtL^-^-prp-^ 


_^ , \ i? 

4n„lp} + P^d) log ij) 

{uj ,/p} 


(3.16) 


where we used uq = 1 and qo = R, and multiplied by (—1)^ , see (1.1). 

It can be shown that for the two special cases p = 0,1, the boundary term above is 
the only contribution to Thus the p = 1 case gives: 


Q{univ) _ 

*^2 — 


(- 1 ) 


d-2 
2 TT 


¥r(i) 


(d+l)r(d-l) 


CtF 


d-3 


(^R^d + qdj 


(3.17) 


^Note that for a CFT, dimensional analysis (and (A.23), (A.26)) dictate the functional dependence of qd 
and Ud such that: qd = and Ud = f 2 {l, n, R, d), where /i, /2 are some functions and tF = nf. 

Note that there is no dependence on mi,... ,mp-\, since the bulk EOM’s (A.23), (A.26) do not depend on 
them. Likewise it can be seen that the sphere result (1.7) doesn’t depend on the m’s. 
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We see from (3.17) that the sign of depends solely on the coefficients and 

of the log term in the FG expansions. These coefficients can be obtained by solving the 
bulk EOMs (A.23) and (A.26). It might be interesting to understand if generally these 
coefficients are constrained to have a definite sign. The p = 0 case (the plane) is examined 
in the following subsection. 

3.2 Plane Entangling Snrface 

The plane entangling surface is a special case of the cylinder with p = 0. Therefore 
(3.16) becomes®: 


S 


(univ) 

2 


, . d — 2 d-|-2 / rJ \ 

(-l)^vr^r(^) 

{d+l)T{d-l) 


{Uj} 


(3.18) 


So the sign of depends solely on the coefficient Ud, which we shall now determine. 

Eor a plane entangling surface and Einstein gravity in the bulk, we will find an explicit 
solution® (A. 19): 


"{rii} 


(z) = jjZ2Kd{hz) 


(3.19) 


where M is a normalization constant, and n 
function gives (see (A.16)) Ud = ■ 


^ small z expansion of this 

Therefore we have: 


Q(univ) _ 

Dr, — 


d-2 
7T 2 


(d-1) 


2'^-2r(d + 2)(f - 1)! 


CtL' 


d-2 


E 

{rij} 




(3.20) 


Since the above expression is positive, we have proved that a plane is a local minimum 
for Einstein gravity in the bulk. It is now natural to conjecture that (3.20) holds for any 
GET. (3.20) was derived for d = even, but it would be simple task to obtain the analogous 
d = odd result. The d = odd result will differ from (3.20) only by its d dependence, and 
for d = 3 it was obtained in [77]. 

Eor the case d = 4, (3.20) gives: 

3 r 2 

n2=l n3 = l 

(3.21) 

This precisely matches (4.13) which we will obtain in the next section for 4d CETs via 
Solodukhin’s formula. 

®Note that (3.16), (3.17), and (3.18) apply to any asymptotically AdS background. 

®This result was also derived in [52] in the context of ’’entanglement density”. 
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4 The Second Order Correction in Field theory 

4.1 4d CFT: Solodukhin’s Formula 

For a d = even CFT there is a universal log term (1.1): 

and as in (1.4), we can expand in small e: 

oo 

g(univ) S{nniv)^k ( 42 ) 

k=0 

Solodukhin’s formula [65] (which applies to a CFT in 4d) is^^: 

= + i (4.3) 

where E 2 is the Euler density, and the integrals are over the entangling surface S, and 

h = Tri^) - \k^k^ = (A:“, - > 0 (4.4) 

The second fundamental form and extrinsic curvature are defined as: 


k 


a 

liv 



r = Tr(fc“ J = , Tr{k^) = 

(4.5) 


where 7 ^ 1 / is the induced metric on the entangling surface. 

We consider shape perturbations that leave the topology of the entangling surface 
fixed, therefore the change in the Euler term above is zero, and of (4.2) will be 

given by the integral of l 2 - 


Q{univ) _ 

br) — 


C4 

2407r 


d^a^h 


(4.6) 


In (4.4) I 2 is always positive [78, 82] and is zero for the sphere. Therefore the sphere 
locally minimizes I 2 and This argument also works for flat entangling surfaces 

(plane, strips) since these have I 2 = 0 . 

In the following, we calculate for several examples. See also [63, 64]. 


• Example 1: Plane Entangling Surface 


The metric in cartesian coordinates is: 

ds^ = dx^ + + dyf (4.7) 

and C 4 are the a and c anomalies in 4d. Both are normalized such that for a real scalar field their 
value is 1. We will sometimes use Ct = instead of C 4 . 
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Consider a plane at x = 0, where y 2 , 2/3 are coordinates along the surface. Now perturb 
it as follows: 


X = 0 + e/(2/2,2/3) 


(4.8) 


The vector normal to the surface is: 

1 


(1 ) ^fy2'! 


lysj 


where 


The second fundamental form at order O(e^) is: 


fvi = dyj (4.9) 


l.a _ 
— 


0 

ify3fy2y3 + fy2fy2y2) 
(/ys/ysys + fy2fy2y3) 


/y2y3 + fy2fy2y2) 

“^^(/ya/yaya + /ya/yaya 

~^fy2y2 

“^/yaya 

~^fy2y3 

~^/yaya 


(4.10) 


starts at order e since the unperturbed plane is flat. Plugging this in (4.4) gives: 

\/lh = h + O(e^) = y {fy2y2 “ fy^ysf + '^fmys (fy2 + .^ya) + (^-H) 

Note that I 2 starts at order e^, therefore the order e correction vanishes even before inte¬ 
gration over the surface. Now Fourier expand the perturbation: 


/(y2,2/3) = E 0-n2n3 cos(n 2 y 2 ) cos(n 3 y 3 ) (4.12) 

n2,n3=l 

and plug (4.11), (4.12) in (4.6): 

f; f; (4.13) 

^42=1 n3=l 

where we used C 4 = and the integral dycos^{ny) = L, where the width of the 

plane is very large L —>■ 00 . So we got a positive result, and therefore the universal term 
of a plane in 4d is a local minimum. We see that (4.13) precisely matches the result (3.21) 
obtained in holography. 

We can compute higher orders of e in (4.2), and we note that for the plane all odd 
terms vanish: ~ order we get: 

^(umv) ^ ^ ^ ^ ^ 13n^ni + 13n|ni) (4.14) 

n2=l n 3 =l 


• Example 2: Sphere Entangling Surface 
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This case was calculated in [63, 64], and we write their result. 
The perturbed sphere entangling surface is: 


r(6»,(/)) = R[l + ef{0,i 


= R 


1 + C ^ ^ m {OA) 


l,m 


where R is the radius of the sphere. Plugging in (4.6) gives: 

l^m 

This matches the holographic result (1.7) [64]. 

• Example 3: Cylinder Entangling Surface 

The metric in cylindrical coordinates is: 

ds^ = + dy^ + 

Consider a perturbed cylinder entangling surface with radius R: 

r{y,(t)) = R + ef{y,4>) 


We get: 

Vlhl, = ^[2/2 + 3/2 + 2/2^ + 8fU^ - R^fl + 8i?2/2^ _ + 2R^fyy 


Plugging f{y, (p) = Y.n,m ^n,m cos{m(p) cos(ny), 


S. 


(univ) _ ^4 


2407r 


dydcp^h 


tt^L 


2 - 5m2 + 2771"^ + (ni?)2(4m2 - 1) + 2{nRY 


Thus > 0 for m >2 and for all n, and the cylinder is a local minimum. 


Example 4: 4d Surface of Revolution 


Consider a surface of revolution entangling surface (e.g. Fig. 2-Left) given by: 

r((9,(/) = ro((9) 

where we use spherical coordinates (r, 0,/)). This surface has rotational symmetry 
(j) direction. Now perturb the surface as follows: 


(4.15) 

(4.16) 

(4.17) 

(4.18) 

(4.19) 

(4.20) 

(4.21) 
in the 


r{9,(t)) = ro{9)[l + ef 2 { 0 ) ■ h{4')] 


(4.22) 










with / 3 ((/>) = x; a-m cos{m(j)). (4.3) and (4.4) give at most 4 derivatives of (j), thus 
the result can be written as a polynomial in m (after integrating over cp): 


tt^Ct 

80 




de 

m^O 


Gi(ro, f2,0)m^ + G2{ro, f2,9)m^ + G^iro, / 2 , 9) 


(4.23) 


where the Gi{ro, f 2 , 9) are some functions of vq, /2 and their derivatives. Explicit calculation 
gives the coefficient: 


Gi{roj2,9) 


tt/I (^4 + 4csc^(0)rQrQ® + 16 csc^( 0 )rQrQ® 


+ 24 csc^(0)rQr(j^ 



8rQ sin0(rQ + 


(4.24) 


Since rg > 0, we see that Gi(ro, f2,9) > 0. Therefore for perturbations with large enough 
m (short wave-length perturbations) (4.23) is positive, and thus all 4d surfaces of revolution 
are local minima. We are not able to show that the functions G 2 (ro, f2,9), G 3 (ro, / 2 , 9) are 
positive (though the final integrated result may still turn out to be positive). 


• Example 5: 4d Waveguide Surface 


Consider a general waveguide entangling surface (e.g. Fig. 2-Right) given by: 

r{yA) = ro{(t)) (4.25) 

where we use cylindrical coordinates {r,y,4>). This surface has translational symmetry in 
the y direction. Now we perturb the surface as follows: 


r = ro(()))[l -h e/2(</')/3(y)] 


(4.26) 


with f 3 {y) = Y, ^^0 cos(my). The result can be written as a polynomial in m (after 

integrating over y): 


tt^Gt 

80 


m^O 


Gi(ro, f2)m^ + G2(ro, f2)w? + G^irg, 72 ) 


(4.27) 


where the Gi{rg, 72 ) are some functions of ro, 72 and their derivatives. Explicit calculation 
gives the coefficient: 


Gi(ro,72) 




2rl^ 


+ 6r§r(? Gr^r'n 2rnr: 


0'0 


6 ^'4 


0'0 


r.4^'6 


'O'O 


4( 


+ r'2 


> 0 


(4.28) 
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Therefore for perturbations with large enough m (short wave-length perturbations) (4.27) 
is positive, and thus all 4d waveguide surfaces are local minima. We are not able to show 
that the functions G' 2 (ro,/ 2 ), G' 3 (?’ 0 )/ 2 ) are positive (though the final integrated result 
g<^niv) positive). 

4.2 Renyi Entropy 

There is a generalization of (4.3) to Renyi entropy in a 4d a CFT [79], [70]: 

4"””’ = <3.(-i)S) + 

(4.29) 

where jg universal log term for the g'-th Renyi entropy. Qi, 2 {q) are functions of 

q such that Qi(l) = < 52 ( 1 ) = 1, so it matches (4.3) in the EE limit q —)• 1. Up to these 
functions, (4.29) and (4.3) are the same, hence the results in the previous sections can be 
used. In particular, if Q 2 {q) is positive then for any entangling surface the Renyi entropy 
will have the same sign as the EE. Note that Q 2 {q) is positive for free helds, and there is 
strong evidence that it is positive for holographic CFTs [70, 82, 89]. If this turns out to be 
correct, then the Renyi entropy will be a local minimum whenever the EE is. 

4.3 Free Massive Fields 

Certain universal EE terms for free massive field theories have been found [16-19]. The so 
called ’’universal area law” for free scalars or fermions with mass m has the following form: 

(—, d = odd 1 

, d = even\ 

where 7 ^ is a positive constant that depends only on the dimension d. All of the shape 
dependence is contained in As, the area of the entangling surface. Table 1 lists the sign 
factors (—1)^ in (4.30) for a Dirac fermion, a conformally coupled scalar, and a minimally 
coupled scalar. 


S = 


(4.30) 



Dirac Fermion 

Conformal scalar 

Minimal scalar 

d = odd 

/ d —1 

(- 1 )^ 

(- 1 ) 2 

/ \ d— 1 

(- 1 )^ 

d = even 

(- 1)2 

(- 1 )— 

(- 1)2 


Table 1. The sign factor in the ’’universal area law” (4.30). 

As an example, let us now compute the ’’universal area law” term for a deformed 
sphere entangling surface: 


r{^d-2) = R 1 f X/ ^{lm}^{lm}{^d-2] 

{Im} 

The area of the deformed sphere is: 


(4.31) 


As = / dD,_2/- = ^ + 0(^3) 


F(^) 


{Im} 


(4.32) 
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where we plugged (4.31) and performed the integrals. 

The first term on the RHS is the area of the undeformed sphere, the 0(e) correction 
vanishes as expected, and the O(e^) correction is positive. Therefore the O(e^) correction 
in S has the same sign as the zeroth order, which can be read from Table 1. Generalizing to 
non-spheres, it is easy to see that the area of a perturbed surface of revolution is larger than 
that of the unperturbed surface of revolution. It would be interesting to perform a similar 
analysis to higher curvature terms for free massive fields (i.e to curvature corrections to 
the ” universal area law”). 

A similar analysis applies to ’’universal area law” terms in interacting theories [19, 
21, 69], with a shape dependence that comes only from the area The EE for a GET 
perturbed by a relevent operator of dimension A = contains the following term: 

d+2 

where A is the coupling constant. Such log terms occur both in odd and even dimensions. 

5 Discussion 

In this work we studied the shape dependence of entanglement entropy by deforming sym¬ 
metric entangling surfaces. We showed that entangling surfaces with a rotational or trans¬ 
lational symmetry locally extremize the EE with respect to shape deformations that break 
some of the symmetry. This result applies to EE and Renyi entropy for any QFT in any 
dimension. Using Solodukhin’s formula and holography, we calculated the 2nd order cor¬ 
rection to the EE for CETs and simple symmetric entangling surfaces. In all cases we found 
that the 2nd order correction is positive, and thus the corresponding symmetric entangling 
surface is a local minimum. Perhaps this result holds more generally for any symmetric 
entangling surface. 

Let us mention some possible future directions. 

• The calculation in section 3 considered only Einstein gravity in the bulk, and it would 
be interesting to consider also higher derivative gravity. For spheres, [64] found that 
S^niv) (depends just on Ct and not on t 2 or (these are the three parameters in the 
3-point function of stress tensors). It would be nice to check if this continues to hold 
for other entangling surfaces. 

• Gomputing the EG coefficients Ud and in (3.16) by solving the bulk EOMs for 
the cylinder (A.23) and (A.26). Maybe it is also possible to show, using a more 
general principle, that qa and Ud in (3.17) must have a definite sign. It would also 
be interesting to consider (in holography) more general entangling surfaces with a 
symmetry. 

• The work of [69] attempted to compute for a plane entangling surface in a d = 4 
GET using the perturbative formalism of EE [66]. They were not able to obtain the I 2 
term of Solodukhin’s formula (4.3), and the current situation is somewhat puzzling. 
Our result (3.20) as well as (1.7) (obtained in [64]) might help in resolving this puzzle. 
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• It would be interesting to repeat the analysis of section 4.1 for 6d CFTs using the 
results of [80, 81]. 

• Other possible extensions are to compute higher orders for j > 2, and also to 

perform computations in a curved space-time background. 
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A Equations of Motion for the Minimal Bulk Surface 

In this section, we obtain EOMs for the bulk minimal surfaces of cylinder, strip, and plane. 
We solve the 1®* EOM for the plane, and obtain Ud in (3.18). 

A.l Strip Entangling Surface 

Consider a strip entangling surface and a bulk metric in Poincare coordinates: 

^AdS r 1 


ds'^ = 


-dz^ — dt^ + dx^ + dyf 


(A.l) 


V(5{z) 

Where I3{z) is some function of the holographic coordinate z, which for AdS: j3{z) = 1. x 
is the direction perpendicular to the strip entangling surface, yi are directions along the 
strip entangling surface, and i = 1,... ,d — 2. The holographic EE is (see (3.1), (3.2)): 

1 /* rzmt 




dz- 




f5{z) L 


1 (dy.x)^ + {dzxY 


(A.2) 




In the ansatz of (3.9) this corresponds to p = 0, Fi = -, F 2 = ld{z), and 

4Gjv\//3(2) 


F 3 = 1. The corresponding EOM is: 


1 


d 


dz \ z'^~^ 


dzX 


^ ^ + {dyiX? +{dzxf' 


+ 


^_^ 

dyi 


dy,x 


^ l + (9j;,x)2 Fidzxf' 


The bulk surface can be expanded in e: 

x{z,yi) = xq{z) + exi{z,yi) -h €^X 2 {z,yi) -t- .. 


= 0 


(A.3) 


(A.4) 


where because of the translational symmetry of the strip, xo{z) does not depend on yi. 
The EOM (A.3) at 0*^ order in e is: 


d 


dzXQ 




= 0 


The solution to this is: 


{d^xof = 


z 


2d-2 


ld{z){Zmax - 


(A.5) 


(A. 6 ) 
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where Zmax is the turning point of the bulk surface. The EOM at 1^* order in e is: 


d / 1 






3/2 


+ 


dlxi 


yd—1 


.P(z) 

Plugging (A.6) in (A. 7): 




= 0 


dz I 

Simplifying, we get: 


^2d-2 

^max 


dlxi + 


^'{z) jd-l) 

2/3(2;) 2; 


1 + 2 


2d-2 


1 - 


2d-2 


dzXl + 


1 - 


2d-2 


(A.7) 


<< I M ^ _ „ (a.8) 


^ = 0 (A.9) 


This equation seems hard to solve analytically. In the next section we consider the 
simpler case of a plane entangling surface. 


A.2 Plane Entangling Surface 

For a plane entangling surface situated at x = 0, the bulk minimal surface goes straight 
down in the bulk: xo(2;) = 0. We have (see (A. 4)): 

x{z,yi) = €Xi{z,yi) + e^X 2 {z,yi) + ... (A.IO) 


The turning point of the bulk surface is at Zmax —^ 00 , therefore (A.9) becomes: 


d^xi + 


/3'{z) jd-l) 

2/3(2;) z 


dzXi + dy.xi = 0 


We Fourier expand xi: 


d-2 


xiiz,y) = ^a{„jx||(.^( 2 ;) JJcos(niyi) 

{ui} i=l 


(A.ll) 


(A.12) 


where the rii are integers, and we defined the shorthand notation: {n*} = ni,... ,nd- 2 - 
Plugging this in (A.ll) gives: 


^ {rii} ^ 


^'{z) [d-l) 

2/3(2;) 2; 


dzX^P 1 — n^x 1=0 

^ {rii} {Ui} 


(A.13) 


where we defined Now we consider the ansatz /3{z) = 1 + az^ for the bulk 

metric, then: 


5^x1^^ 1 + 

^ {rii} 


akz^ 


-1 


d-l 


2(1 + az^) z 


a ( 1 ) -2 ( 1 ) 

OzX) 1 — n Xr 1 
^ {rii} {rii} 


= 0 


(A.14) 


- 17 - 




























(A.15) 


The CFT case j3{z) = 1 can be recovered by plugging a = 0: 

dl 


(1) _ ^ (1) _ ^2 (1) _ Q 


This equation is solved by modified Bessel functions^^: 


\{^) — Ciz‘2 K d(hz) + C 2 Z 2 Id(hz) (A.18) 

2 2 

Choosing a boundary condition such that does not explode at z —)■ 00 , leaves 

only Kd. Therefore the final solution is: 

2 

where the normalization constant is M = 22“^(| — l)!n “2 . We will use the above solution 
in (3.19). This result was also derived in [52] in the context of ’’entanglement density”. 


A.3 Cylinder Entangling Surface 

Consider a cylinder entangling surface x in flat space-time M'’*. The holographic 

EE in cylindrical coordinates is (see (3. 8), (3. 9)): 

^ ^ / dzd^-^-^ydn,^^^jl + (5,r)2 + {dy^rY + ^(dn,r)2 (A.20) 

where we assumed the AdSd+i metric. The corresponding EOM is: 


d 


dz \ z‘^~^ 


dzV 


+ 


1 d 

yd—1 




1 d 


yjl + (dzry + (dy.r)^ + ^{dn^rY) ^ V \Jl + (9^r)2 {dy.rY + ^{dn^rY 

\ 1 d 


rP 


dn. 


'p V Y^l + {dzvy + {dy^Y + 


zd-i dr I 1 + ;T(^f2,r)2 ] = 0 


(A.21) 


The bulk surface expanded around the 0*^ order cylinder is: 

r = ro{z) + eri{z,yj,Qp) + d‘r2{z,yj,VLp) + ... 

The modified Bessel function: 

KAz) = ^ {lY'' E(-ir ~ + {-ir+Aog{z/2)Uz) 

I / 1 ^0' + 1 ) + ^(^ + i + 1 ) ...2J 

^ ^ 2 I 2 J ^ Tj!(A+j)!^ 

where ip is the digamma function, and 


(A.22) 


(A.16) 


(A.17) 
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The 0*^ order EOM is: 


d ( rl dzTo 


dz ^ ) z 

The 1®* order EOM is: 


- -i[=iP^l Vl + (5^^o)^ = 0 


(A.23) 


d^ri + 


rAi To 


dz\yz‘i ^ (l+(9zro)^)^T 


z'^ 1 {l+{dzroY)^/^ 


dzVi + [l + {dzrof] (+ ajri) = 0 


(A.24) 


We Eourier expand ri (see (3.11)): 


d—2—p 


ri{z,yj,np)= n cos(niyi) (A.25) 

{'Tlj ^Ip } ^ 1 

where the rij are integers, and we defined the shorthand notation: {nj} = ni,..., n^_ 2 _p 


and {Ip} = I, mi ... We plug (A.25) in (A.24): 


j 1 + 

2 {nj,lp} 


rAlTo 


dz yz'^ 1 (l+(9zro)^)®T 


z‘‘ 1 (l+(9zro)^)^T 


+ [' + 


p-i(,+p-l) _,,N (1) 


{nj,lp} 


(A. 


where we defined nj, and used d'^Yi^ = —l{l+p — l)Yi^ [64]. Since the parameters 

mi,... ,mp_i do not appear in (A.26), the solution ^ | will not depend on them (and 
neither will *§’ 2 °“"''^'|in (3.16)). 

Eor a plane entangling surface we have p = 0 and r = x, and d^XQ = 0, and in this 
case (A.26) reduces to (A. 15). On the other hand, for a sphere entangling surface we have 
p = d — 2 and Tq = B? — z^. Thus for a sphere the EOM is: 


9.2 


,(i) 


z' l 


l {d-l)R^ + 2z^ ( 1 ) R^[d-2-l{l + d-^)] ji) ^ 

Z R^-Z^ " ' (i ?2 _ ^ 2)2 


(A.27) 


We find a solution to this equation in terms of hypergeometric functions. Eor d = 3 
the solution is: 


( 1 ) ^ 1 / z- R \ 2 / R + lz \ 

' M\z + r) 


(A.28) 


which agrees with Eq. 43 of [63]. 
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